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ABSTRACT 

The concentration-mass relations proposed by Prada et al. (2012) and by Duffy et al. 
(2008) on the scales of galaxy clusters show some of the largest discrepancies among 
all the works present in literature. This is surprising because they are both derived 
from the analysis of dark-matter halos forming in ACDM simulations with similar set- 
ups. With the help of analytic models and numerical simulations we investigate the 
origin of this discrepancy focusing on the procedures used to derive the concentrations 
(circular velocity ratios versus density profile fitting) and on the selection criteria 
used to bin the halos (binning in maximum circular velocity versus binning in mass). 
We find that both steps of the analysis have a large impact on the resulting c — M 
relation. In particular, we show that the two c— M relations can be entirely reconciled, 
if we account for these methodological differences. Our analysis demonstrates that 
the concentration estimates are sensitive to the largely different radial scales probed 
by a particular measurement method. This implies that concentrations derived with 
different techniques (both in observations and in simulations) must be compared over 
the same radial range. 

Key words: galaxies: clusters - cosmology: theory - dark matter - methods: analyt- 
ical, numerical 



1 INTRODUCTION 



About fifteen years ago Navarro et al. ( 1996 1 showed that 



dark-matter halos forming in simulations of hierarchically 
clustering universes develop a nearly universal density pro- 
file, reasonably well described by the formula 

Ps 

R/R 3 (l + R/R s ) ' 

where p s and R a are the halo characteristic density and 
scaling radius, respectively. This result holds over a wide 
range of masses, from the scale of massive galaxies to that 



p(R) = 



(1) 



of galaxy clusters. Navarro et al. ( 1996 1 also found that 



low mass halos are centrally denser than their high-mass 
counter-parts. In particular, the halo concentration, defined 
as C200 = R200/ Rs, was found to be a decreasing function of 
the halo mass, M2ocQ This outcome was interpreted as the 
signature that a dark matter (DM) halo keeps memory of 
the density of the universe at the epoch of its collapse. For 
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1 ^200 an d M200 the radius and the mass of the sphere en- 
closing a mean density equal to 200 times the critical density of 
the universe. 



this reason, the c — M relation is connected to and provides 
information on the cosmological model within which halos 
form ( jNavarro et al.|1997 |. 



Virtually, all cosmological observations of the universe 
on large scales undertaken in the last decade support the so- 
called concordance cosmological model (ACDM). Simulated 
universes reproducing variations of this model have been 
performed to calibrate the c — M relation (Navarro et al.| 



tacharya et al. 



2011 



T9971 |Bullock et al.|[200l| |Eke et al.|2001| [Zhao et al.||2003 



Duffy et al.|2008||Gao et al.|2008||Macci6 et al.|2008|fBhaTl 



Klypin et al. 20111. Most theoretical 



works show relatively similar behaviors with the exception 
of the recent analysis of Prada et al. (20121, which differs 



form the others both in terms of amplitude and shape of the 
c — M relation, especially on the scales of galaxy clusters. 

In this paper, we investigate the origin of this discrep- 
ancy focusing on the comparison between the relations of 



Prada et al.| ( [2012| P12) and |Duffy et al.| ( |2008| D08). Both 
works are based on the analysis of a large number of DM 
halos extracted from cosmological boxes. For the scope of 
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studying clusters, P12 recurred to the Multidark]^] simula- 
tion consisting of a l/i" 1 Gpc 3 box filled with 2048 3 DM 
particles. Instead, D08 used smaller boxes, the largest one 
having a side-length of 400/i _1 Mpc filled with 512 3 par- 
ticles. These simulations are therefore characterized by a 
worse resolution and smaller statistical power on the clus- 
ter scales. The simulation codes have different architectures 
(mesh-based, ART 



1997 



et a 



in the first case 



|2001| in the sec- 



Kravtsov et al. 
and particle-based, GADGET - |Springe 

ond). However, since a long time ago, it has been demon- 
strated that the DM properties are unaffected by the nu- 
merical scheme ( Frenk et a.1. || 1999 > . Finally, the differences 
in terms of cosmological parameters are smalQ 

That said, it is surprising that the P12 concentrations 
are ~ 50 — 60% higher than the concentrations found by D08 
at fixed halo mass. The P12 relation, additionally, exhibits 
an upturn towards the largest masses: after a certain "pivot" 
mass the logarithmic slope of the relation from negative be- 
comes positive. Such behavior is particularly evident at high 
redshift, as the pivot mass decreases with redshift. 

Some discussion on the impact of selection effects on the 
shape of the P12 c — M relation can be found in other works 



(e.g. Ludlow et al. 20121. Here, we aim at fully shedding 
light on the reasons of the mismatches, discussing the dif- 
ferences between the two methods employed by P12 and by 
D08 to measure concentrations and to construct the c — M 
relation (Sect. [2]). To be immune from numerical nuisances, 
in Sect. [3] we study the systematics introduced by the two 
procedures utilizing both analytic models and the halos from 
the best resolved Multidark simulation. Finally, in Sect. [4] 
we summarize the results and draw our conclusions. 



at R„ 



Duffy et al. 2008, instead, followed a more standard 
methodology. They searched for the NFW best-fit of the halo 
density profiles (Eq. 1) evaluated in equally logarithmically- 
spaced radial bins, with —1.25 ^ log R/ R V i r ^ 1. The con- 
centration (cdos) is one of free parameters of the fit (the 
other is the normalization p s , Eq. [T]) . 



2.2 Binning 

On top of the concentration derivation, the procedures of 
P12 and D08 differ also for the selection criteria adopted to 
bin the halos. In the former case, halos are sorted by Vmax 
(cvsei), while in the latter via mass (cmsci)- 

This has important consequences. Using Eq. [2] and the 
circular velocity definition, we can see that 



Vmax = V 200 J°-^ = ] f^F(c) oc M^F(c) , (3) 



Thus, binning in Vmax translates into selecting curved re- 
gions in the concentration-mass plane as shown in Fig. 1 
where some curves at constant Vmax (from 1000 km/s to 
1500 km/s) are plotted. Note that the curves are steeper 
at the high-mass end while they flatten towards the smaller 
masses. Consequently, we expect the median concentrations 
in V max bins to be larger than those measured in mass bins, 
because the same V max bin also contains many systems hav- 
ing slightly smaller mass but much higher concentrations. 



2 PROCEDURAL DIFFERENCES 

2.1 Methods to derive the concentration 

Prada et al. (2012) 's procedure is described by the fol- 
lowing steps: 

1) each halo is identified in the cosmological box and char- 
acterized by means of its -R200 radius and M200 mass; 

2) two velocities are computed from the mass profile in log- 
arithmically equi-spaced radial bins: 

2) the halo maximum circular velocity: V max = 
max(V C i r c), where V c i rc = \/GM{< R)/R, and 

ii) the circular velocity at -R200: V200 = y GM200/ ' R-20Q- 
Assuming a NFW profile, the two velocities are related 
through the concentration in a way that more centrally con- 
centrated halos are characterized by higher V ma x /V200 ra- 
tios: 



V„, 



V20O 



0.216c 
~1W 



= F(c) , 



(2) 



where /(c) is given by /(c) = ln(l + c) — c/(l + c); 
3) finally, the halo concentration is derived by numerically 
inverting Eq. [5] P12 concentration (cpia) is, therefore, 
estimated from mass measurements done only at R200 and 



2 http://www.multidark.org 

3 (n m ,Cl b ,n A ,h,a B ,n a ) equal to (0.27, 0.046, 0.73, 0.70, 0.82, 
0.95) for P12 and (0.26, 0.046, 0.74, 0.80, 0.80, 096) for D08 



3 UNDERSTANDING THE DIFFERENCES 
3.1 Consequences of the binning 

To illustrate the effect of the binning method, we use ana- 
lytic prescriptions and Monte-Carlo simulations to generate 
halos drawing them from the Sheth & Tormen mass function 
| |Sheth fc~T ormcn 2002) and associating them the fiducial 
concentration c(M, z) of D08 for relaxed clusters: 



C200(M200, z) 



6.71 



Ah 



(1. + z) 0A4 \2x lO^h^M® 



(4) 



Further, we assume that the concentrations at fixed mass 
are log-normally distributed, 



p(c) = 



: exp 



(lnc-lnc) 2 



2a 2 



(5) 



and we change the concentration scatter by varying a 
between 0.2 and 0.6 in steps of 0.1. Numerical simulations 
are usually consistent with a ~ 0.25 — 0.4 ( jDolag et al.| 



2004 De Boni et al. 20121, almost independently of mass. 



We find that the Multidark halos used by P12 have a 
concentration scatter of a ~ 0.35. To focus exclusively 
on the problematics of the binning, we do not measure 
the concentrations of these halos a la P12, but we adopt 
their binning method to re-construct the c — M relation. 
We verified that the choice of another mass function or 
concentration-mass relation to generate the halos does not 
change the results of the following analysis. 
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Halos at low redshift. According to previous equa- 
tions, we begin our analysis by building concentration dis- 
tributions for different values of a and for halos at redshift 
z = 0.25. First, we select halos in a narrow bin of circu- 
lar velocity: 1000 km/s < V m „ ^ 1100 km/s. The median 
mass in the bin varies little as a function of the concen- 
tration scatter, being 2.1 x 10 14 h~ 1 M t r ) for a — 0.2 and 
1.9 x 10 14 ft _1 M for a = 0.6. Using Eq.|iJ these mass val- 
ues correspond to D08 concentrations of CMsei ~ 3.96 and 
CMsei ~ 3.99 respectively, but due to the positive skewness of 
the concentration distributions in the Vmax bin, we measure 
median concentrations for the selected halos growing from 
CD0$,Vsei = 4.0 for a = 0.2 to C£>o8,v,»ei =5.2 for a = 0.6. 

These results confirm that median concentrations are 
larger in Vma^-selected than in mass selected samples of ha- 
los, even if the median mass in the bin remains (almost) con- 
stant. Moreover, the ratio between median concentrations in 
Vmax- and mass-selected samples is a growing function of the 
concentration scatter of the parent halo population. 

In the left panel of Fig. [T] we generalize these results 
for the entire population of halos assuming a concentration 
scatter a = 0.4. The medians of the concentrations in each 
Vmax bin are shown as black diamonds with the associated 
inter-quartile ranges (vertical bars). The blue dashed line 
represents the best-fit to the data using the P12's functional: 



C200 = a 



M 2 qq 
lO^/i-iMo 



1 + 7 



/ M200 
VlO^h-iMg 



1/2 



(6) 

The original relations by P12 and D08 are shown by the 
dotted-magenta and by the solid-red lines respectively. 
Finally, as a consistency check, the concentrations derived 
by binning the sample in mass are plotted with red crosses. 
As expected these data points perfectly reproduce the 
input D08 c — M relation. As discussed above, we find 
that for any mass the concentrations measured in Vmax 
bins are higher than those derived in mass bins. The ratio 
of concentrations in Vmax (cvsel) an< A in mass bins (cMsei) 
grows with the concentration scatter as shown in the central 
panel of Fig. [I] It also depends on the halo mass, being a 
slowly decreasing function of M200 at low masses, but with 
a clear upturn at the highest masses. For a ~ 0.3 — 0.4, 
typical of numerically simulated halos, the concentrations 
of the VWiai-selected halos are higher than those of the 
mass-selected ones by an amount that grows from 5% to 
20% between M = 10 14 /i _1 M and M = 10 15 /i _1 M Q . At 
the largest masses, the ratio is strongly dependent on a, 
and it changes from ~ 5% to ~ 55% as a grows from 0.2 to 
0.6. 



the effects of binning. Results for lower values of a can be 
extrapolated from the previous discussion. 

We clearly see an upturn in the c—M relation when 
halos are selected by Vmax- The upturn begins at smaller 
masses as the redshift increases. Interestingly, we see that 
binning in Vmax also inverts the redshift evolution of the 
input D08 c—M relation (solid lines). On the cluster scales, 
the concentrations of V ma a:-selected halos grow with redshift 
at fixed mass, reproducing the trend which characterizes the 
P12 c—M relation. For 2 = 0.5 and a — 0.6 binning by 
Vmax is sufficient to bring the D08 and the P12 to a full 
agreement. However, we notice that it becomes harder and 
harder to reproduce the P12 relation as the redshift increases 
even assuming an extreme concentration scatter. In order to 
obtain a good match with P12 at 2 > 0.5, we would need 
to modify the initial c — M relation of D08 by increasing its 
amplitude and reducing its redshift evolution. We verified 
that this can achieved by changing the exponent of the (I+2) 
term in Eq. [4] from 0.44 to 0.25. 

In the literature, there is no consensus on the evolution 
of the c—M relation with redshift, in particular for massive 
halos. While several authors predict a strong redshift evolu- 



tion of the concentrations at all mass scales (see e.g. Bullock 
|et al.|2001||Eke et al.|200"Tj ), |Zhao et~aT1p003| find that the 
evolution of the concentration of individual halos is not just 
a function of redshift but is tightly connected to their mass 
growth rate (see also |Wechsler et al.| [2002). In particular, 
the faster the mass grows, the slower the concentration in- 
creases. Since most of the massive halos are in a fast mass 
accretion phase at high redshift, they find that the cluster 
c—M relation has a very slow redshift evolution. Similar 



results were found recently by Munoz-Cuartas et al. (2011 1, 



who confirm that the growth rate of the concentration de- 
pends on the halo mass, with low-mass haloes experiencing a 
faster concentration evolution. These authors also find that 
the evolution of the c — M relation is faster at lower redshifts 
than at higher redshifts. 

To summarize this section, we showed that binning in 
Vmax instead of binning in mass has two consequences: 1) 
the median concentrations are higher, and 2) there is an up- 
turn at the large mass scales. Both these effects are amplified 
by increasing the scatter in the concentration distribution. 
Nevertheless, the difference in the binning procedure cannot 
be the only cause of the discrepancy between P12 and D08. 
Indeed, as shown in Fig.[T] for values of a that are typically 
found in the simulations the increment in the concentration 
value caused by binning in Vmax is ~ 15 — 20%, which ac- 
counts of about half of the differences between the P12 and 
D08 c-M relations. 



Halos at higher redshift. P12 showed that the up- 
turn in their c — M relation becomes more prominent at high 
redshifts. We repeat the same procedure to generate halos 
at redshifts 2 = 0.5, 1,2 and 3 (although the D08 relation 
has been calibrated up to 2 = 2). The results in these cases 
are shown in the right panel of Fig. [T] To avoid confusion in 
the figure, we do not show the data points here, but we plot 
only the best-fit curves, given by the dashed lines, whose 
thickness increases with redshift. For this analysis, we con- 
sider only an extreme concentration scatter, a — 0.6. We do 
not claim that this is realistic, but we use it to emphasize 



3.2 Consequences of the concentration 
measurement 

Now, we investigate if the remaining discrepancy can be at- 
tributed to the different methods employed by P12 and D08 
to measure the concentrations (Sect |2.f| . To do this, we use 
numerically simulated halos and, to avoid any extra bias 
introduced by using different samples, we analyze halos ex- 
tracted from the Multidark simulation, the same simulation 
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Figure 1. Left panel: The c — M relation of V mtl:E (diamonds) and mass-selected halos (crosses) at 2 = 0.25. Results are shown for 
cr = 0.4. The solid-red and the dotted-blue lines represent the D08 and P12 c — M relations, respectively. The dashed-blue line indicates 
the best-fit to the diamond data points. We also show some curves of constant V max (thin solid lines). Central panel: For different values 
of cr, ratio between the concentrations of halos selected by either Vmax or mass as function of mass. Right panel: Evolution of the c — M 
relation of V maa: -selected halos as a function of redshift (dashed lines). Results are shown for 2 = 0.5, 1, 2, 3 and refer to the case a = 0.6. 
Ticker lines (red-yellowish colors) correspond to higher redshifts. The solid and the dotted lines show the redshift evolution of the D08 
and of the P12 c — M relations. 



used by P12. From their databas^J we query for the masses, 
concentrations (measured from the V- ratio), maximum cir- 
cular velocities, and mass profiles of all halos at redshift 
z — 0.25 with masses satisfying at M200 ^ 10 14 /i _1 Mq. 

To corroborate the final concept of the previous section, 
we create a plot similar to the left panel of Fig. [I] where the 
concentration are now derived following the method of P12 
(Fig. [2J. Black diamonds and red crosses still denote re- 
sults from the binning in V ma x and mass, respectively, while 
dotted-magenta and dashed-red lines report the original re- 
lations of P12 and D08. As consistency check, we can notice 
that the P12 relation is perfectly recovered. Selecting ha- 
los by mass, we obtain concentrations lower by ~ 5 — 15% 
compared to P12. However, this is still ~ 15 — 20% higher 
than the amplitude of the D08 relation. Once again, even 
if we change the method to derive the concentration, the 
total discrepancy between D08 and P12 is not explained by 
simply modifying the binning criterium. As final step, from 
the mass profile of the Multidark halo we derive the concen- 
tration a la D08 and plot the results after binning in mass 
(orange crosses). These values are now consistent, within ler 
with the original D08 relations shown by the yellow band 
whose upper and bottom limits refer to relaxed clusters and 
un-relaxed halos. The fact that they better match the upper 
limit may be due to the better resolution and larger size of 
the P12 simulations. 

The circle is now closed, we demonstrated that the two 
steps in the procedure to construct the c — M relation are 
concurring to produce a difference of ~ 40% between the 
results of P12 and D08. This result being established, we still 
need to answer to the question: "why do different methods 
to derive the concentration lead to discrepant results?" 

To discuss this issue we extend the experiment pre- 
sented in P12 where they compare the concentrations mea- 
sured from the V = Vmax/Vzoo ratio and from the density 
fit of the 10000 more massive systems in the Bolshoi simu- 
lation at z — and at 2 = 3. P12 found that the agreement 
between the concentrations is good for the highest values 
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Figure 2. The c — M relations of cluster-sized halos extracted 
from the Multidark simulation at 2 = 0.25. The black diamonds 
show the result of selecting clusters by Vmax- The P12 c — M 
relation is indicated by the dotted line. Red crosses show the re- 
sults obtained by binning halos in mass. The best-fit relation is 
given by the dashed line. Finally, the orange crosses show the 
median concentrations in mass bins, when the concentrations are 
measured by fitting the mass profiles with NFW models, i.e. fol- 
lowing the approach of DOS. The shaded yellow region shows the 
D08 c — M relation (relaxed and un-relaxed halos). 



of c, i.e. for the least massive halos. For c < 5, instead, the 
concentrations cpi2 are higher than the cdos by ~ 10 — 15%. 

Using the sample of halos with M > 10 14 /i _1 Mq at 
z = 0.25, we find that, if we fit the density profiles with the 
NFW functional in Eq. [I] over the radial range [R m ax , R200] , 
i.e. in the same radial range probed by the concentrations 
cv of P12, the resulting concentrations c/« are on average in 
good agreement with those derived from Eq. [2] The results 
are shown in the left panel of Fig. [3] as function of the \ 2 01 
the best NFW fif] 
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5 The fit is done using the IDL routine MPFITFUN - s 
http: // cow. physics .wise . edu/~craigm/idl/mpf ittut .html 



On the cluster concentration-mass relation 5 




10" 6 10" 4 10~* 10° 10* 0.1 1.0 10.0 100.0 10' 4 10' 5 

X 2 X 2 » [h"' M ] 



Figure 3. Left panel: Ratios between the concentrations measured using the V-ratio, as done by P12, and by fitting the halo density 
profiles with the NFW model in the radial range [R ma x — -f?20o]- The ratios are shown as a function of the x 2 value of the best fit NFW 
models. Central panel: as in the left panel, but the fit extends to 0.05R v i r (as done by D08). Right panel: quality of the best fit as a 
function of the halo mass. 



However, we notice that, for an NFW halo, R ma x ~ 2i? s 
I jNavarro et al.|[l997[ ). Thus, the radial range we just used 
is significantly different from the radial range within which 
D08 measured their concentrations. Indeed, they extended 
the fit to much smaller radii (down to 5% of Rvir)- Repeating 
the fit over the radial range used by D08, we found that 
only the halos that are better fitted by an NFW model (i.e. 
having the smallest x 2 ) have Cfl t ~ cv (central panel of 
Fig. [3). This is not surprising, given that Eq. [2] holds only 
for NFW halos. On the contrary, there is a general tendency 
to measure cat < cy for the bulk of halos. The right panel 
of Fig. [3] shows that x 2 is positively correlated with the 
mass. On the basis of these results, we must conclude that 
for a large fraction of halos, and for the most massive in 
particular, the NFW functional does not represent a good 
fit to the density profile. As a consequence, the value of the 
resulting Cfu is dependent on the radial range chosen for the 
fit. 



4 SUMMARY AND CONCLUSIONS 

In this paper, we considered two among the cluster-scale 
c — M relations proposed in the literature that exhibit the 
largest mismatch: those proposed by Prada et al. (2012) and 
Duffy et al. (2008). These relations differ mainly in two 
aspects: i) the amplitude of the D08 relation is lower by 
~ 40%; ii) the c — M relation of P12 exhibits an upturn 
at the highest masses which is not present in the D08 re- 
lation and in other works. By means of both analytic and 
numerical models, we demonstrated that the mismatch can 
be explained by differences in the procedures for binning the 
halos and for measuring the concentrations. In more details, 
we found that: 

• if the c—M relation is constructed using halos binned by 
their maximum circular velocity (P12) higher concentrations 
are expected at the same mass scale than for mass-selected 
halos (D08). This is true independently of the method fol- 
lowed to measure the concentrations. We found that the 
difference between CMsei and cvsel grows with the scatter of 
the concentration distribution. For values of the scatter typi- 
cally found in numerical simulations (a < 0.4), the expected 
impact of the halo selection method on the c — M relation is 
~ 15 — 20% and it is mass-dependent: the highest masses are 
mostly affected by the binning method (see also appendix of 



P12). Binning can also explain, at least partially, the upturn 
seen in the P12 c — M relation; 

• the concentrations derived by fitting the halo density 
profiles with NFW models (Eq. [T] D08) and by means of 
the V^-ratio (Eq. [2] P12) agree well with each other only 
if the NFW fit is performed over the radial range R max ^ 
R ^ -R200, which corresponds to the scales where the V-ratio 
is evaluated. When we fitted the Multidark halos over the 
same radial range used by D08, we found that the agree- 
ment between concentration measurements is good only for 
the halos with the smallest best-fit x 2 ■ However, the bulk 
of halos shows a tendency for larger concentrations when 
using the P12 method. We also found that the best-fit x 2 
is positively correlated with the halo mass. Typically, the 
most massive halos have a worse NFW-fit because they are 
generally far from equilibrium. This helps further to explain 
the upturn in the c — M relation of P12. We verified that the 
method to measure the concentrations impacts for another 
~ 20% on the amplitude of the c—M relation for halos at 
z ~ 0.25. 

On the basis of these results, we can conclude that the re- 
lations proposed by P12 and D08 can be fully reconciled 
if one takes into account the impact of both the halo se- 
lection and the method to measure the concentrations. An 
important caveat that emerges from our analysis is that, to 
compare results, it is relevant to measure concentrations in 
a consistent manner. This is even more crucial when we aim 
at comparing theoretical and observational results. Indeed, 
different observations probe the halo mass profiles over dif- 



ferent radial ranges < |Comerford fc Nataraja n 2007; Ettori 
et al |2011||Schmidt fc Allen|2007| |Oguri et al'|2012| |FedeIi 



2012 1. To this issue we dedicated an entire companion paper 
( Rasia et al.|2013[ ). 
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